Best 

Available 

Copy 


ACCEPTANCE  ^ONTROL  IjgARTS  WITH 
JIROR  PROBABILITIES 

5ASED  on  Bisson  jountj^ata# 


Research  ^ep‘o?t#(^o. 
~~  by 


80-9 


Surcsh^fhatre/ 

'  j  Richard  L./Scheaf ferfr* 

_ _ -/  Richard  S. /Leavenworth^ 


W-1  ! 


*Depnrtment  of  Industrial  and  Systems  Engineering 
University  of  Florida 
Gainesville,  Florida  32611 


**Dcpartment  of  Statistics 
University  of  Florida 
Gainesville,  Florida  32611 


K.rr-WJ-! 


•  ‘  ’  'f-  f 

Am*  t- 


V. 


FEB  i  7  1981 


r,v  v  v 
ilv  Y  •* 

Uj 


A 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


This  research  was  supported  by  the  U.  S.  Department-aL-the  Navy  r 
Office  of  Naval  Research  under  ContrnctT N00^14-75-C-Q'781^y7 

/q 


__ _ w 


THE  FINDINGS  OF  THIS  REPORT  ARE  NOT  TO  BE  CONSTRUED  AS  AN  OFFICIAL 
DEPARTMENT  OF  THE  NAVY  POSITION,  UNLESS  SO  DESIGNATED  BY  OTHER 
AUTHORIZED  DOCUMENTS, 


fCCgmTV  CLASSIFICATION  OF  THIS  FACE  (Whtt  Pit*  Bntm4) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETINO  FORM 


UToST  number 

80-9  ^ 


|2.  JOVT  ACCESSION  NOJ 

i 


m=AM 


*.  RECIPIENT'S  CATALOG  NUMBER 


TITLE  (mi  Submit) 

Acceptance  Control  Charts  with  Stipulated 
Error  Probabilities  Based  on  Poisson  Count  Data 


6Xi 


S.  TYRE  OF  REPORT  S  PERIOD  COVERED 

Technical 


f.  PERFORMING  ORO.  REPORT  NUMEER 

80-9 


»•  CONTRACT  0*  OR  ANT  NUMBERf*) 


N00014-75-C-0783 . 


AUThORTa) 

Suresh  Mhatre 
Richard  L.  Scheaffer 
Richard  S.  Leavenworth 


PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Industrial  and  Systems  Engineering  Department 
University  of  Florida 
Gainesville,  Florida  32611  ^ 


to.  PROORAM  ELEMENT.  PROJECT,  TASK 
AREA  •  WORK  UNIT  NUMBERS 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Office  of  Naval  Research 
Arlington,  VA 


IZ.  REPORT  DATE 

December,  1980 


tv  NUMBER  OF  PAGES 

32 


U,  MONlToAlNO  AGENCY  NAME  S  ADDRESS^  iiUtrtnt  Irom  Controlllni  Olllctj 


IS.  SECURITY  CLASS,  (tl  tfi It  refer t) 

UNCLASSIFIED 


IS*.  OECL ASSI F|C A TIOn/ OOWNORAOINO 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (o I  (lilt  ktpotl) 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED. 


IT.  DISTRIBUTION  STATEMENT  (ol  the  abetted  entered  In  Block  20,  ll  illlerent  tram  Report) 


N/A 


IS  supplementary  notes 


IS.  KEY  WOROS  (Continue  on  rtrtttt  tide  II  noctteary  mi  Identity  by  bloik  number) 

Control  Charts,  Poisson  Distribution,  Normal  Approximation 


20.  ABSTRACT  fCondnu*  on  rtrtttt  tldt  II  ntcttttry  and  Identity  by  block  number) 

An  acceptance  control  charting  scheme  is  investigated  for  the  case  in 
which  observations  consist  of  a  number  of  nonconformances  seen  when  a  process 
is  observed  for  a  certain  fixed  length  of  time.  The  counts  are  assumed  to  have 
a  Poisson  distribution.  Two  normal  approximations  for  finding  the  optimum  sample] 
size  and  control  limit  are  compared  to  the  exact  values  found  through  the  use  of 
Poisson  (or  Chi-square)  probabilities.  Recommendations  for  practical  usage  are 
made  as  a  result  of  a  numerical  study. 


dd 


F?:  n  1473 


EDITION  OF  I  NOV  •*  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  this  PAGE  (IWi»n  Date  Entered) 


1 

,  TABLE  OF  CONTENTS 


PAGE 


ABSTRACT  .  i 

SECTIONS 

INTRODUCTION  .  1 

PROBLEM  FORMULATION  .  2 

PROBLEM  SOLUTION  .  A 


Exact  Solution  .  . . 5 

Standard  Normal  Approximation  ....  .  .  .  5 

Square-Root  Normal  Transformation  .  6 

NUMERICAL  STUDY  .  7 

EXAMPLE  APPLICATIONS  .  10 

Example  1.  Check  of  Repaired  Items  Against  a  Standard . 10 

Example  2.  Establishing  a  Standard  Plan  to  Check 

Maintenance  Errors  in  a  Paint  Shop . 11 

CONCLUSIONS . 16 

REFERENCES . 17 

APPENDIX  -  COMPLETE  RESULTS  OF  NUMERICAL  STUDY  .  18 


ACCEPTANCE  CONTROL  CHARTS  WITH  STIPULATED  ERROR 


PROBABILITIES  BASED  ON  POISSON  COUNT  DATA 

by 

Suresh  Mhatre 
Richard  L.  Scheaffer 

,  Richard  S.  Leavenworth 

I 

i 

i 

*  ABSTRACT 

\ 

\  ' 

An  acceptance  control  charting  scheme  is  investigated  for  the  case  in 
which  observations  consist  of  the  number  of  nonconformances  seen  when  a 
process  is  observed  for  a  certain  fixed  length  of  time.  The  counts  are 
assumed  to  have  a  Poisson  distribution.  Two  normal  approximations  for 
finding  the  optimum  sample  size  and  control  limit  are  compared  to  the  exact 
values  found  through  the  use  of  Poisson  (or  Chi-square)  probabilities.  Rec¬ 
ommendations  for  practical  usage  are  made  as  a  result  of  a  numerica1  study. 


ACCEPTANCE  CONTROL  CHARTS  WITH  STIPULATED  ERROR 
PROBABILITIES  BASED  ON  POISSON  COUNT  DATA 


by 

Suresh  Mhatre 
Richard  L.  Scheaffer 
Richard  S.  Leavenworth 

INTRODUCTION 

The  application  of  control  chart  methods  to  accept  and  reject  the 
output  of  a  process  has  been  described  in  the  literature  on  several  oc¬ 
casions.  Winterhalter  (1945)  suggested  the  use  of  what  he  called  reject 
limits  in  conjunction  with  the  usual  control  limits  to  control  a  process 
average,  X.  So  long  as  process  dispersion  was  held  in  control  and  the 
control  limits  lay  within  the  reject  limits,  virtually  all  product  would 
meet  specifications. 

Hill  (1956)  expanded  this  idea  by  employing  the  reject  limits  in  place 
of  standard  control  limits  for  those  cases  in  which  the  difference  be¬ 
tween  the  upper  and  lower  specification  limits  (U  -  L)  substantially  ex¬ 
ceeded  the  natural  tolerances  of  the  process,  6  o'.  In  his  1957  article, 
Richard  Freund  gave  more  form  and  substance  to  Winterhalter's  earlier 
work  by  providing  an  analytical  basis  for  deriving  the  location  of  a  reject 
limit.  He  also  coined  the  phrase  Acceptance  Control  Chart  and  referred  to 
the  derived  limit  as  the  Acceptance  Control  Limit  (ACL).  His  development 
closely  follows  that  for  variables  acceptance  sampling  plans.  Essen¬ 
tially,  it  requires  the  specification  of  two  points  on  an  operating 
characteristic  curve  in  terms  of  a  quality  level  and  probability  of 
acceptance  for  each.  From  bese  inputs  are  derived  a  control  limit, 
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the  ACL,  and  a  subgroup  size,  n.  So  long  as  the  plotted  values  of  X 
fall  within  the  ACL's,  the  process  may  be  assumed  to  be  turning  out 
product  that  meets  specifications,  subject  to  the  defined  risks. 

In  this  paper,  we  extend  Freund's  work  to  the  attributes  case  of 
counts  of  nonconformances  or  nonconformances  per  unit  which  can  be  shown 
to  follow  the  Poisson  distribution.  Traditionally,  the  £-chart  has  been 
used  when  the  area  of  opportunity  for  a  nonconformity  to  occur  is  constant; 
the  ij-chart  has  been  used  when  the  area  of  opportunity  varies  from  sub¬ 
group  to  subgroup.  Three  methods  for  finding  the  optimum  subgroup  size 
and  acceptance  control  limit  are  compared.  These  are:  (1)  the  exact 
method,  employing  the  Poisson  distribution;  (2)  the  standard  normal  approx¬ 
imation;  and  (3)  the  square-root  normal  transformation. 

i 

PROBLEM  FORMULATION 

Control  charts  for  nonconformances  have  found  many  uses  in  industry. 
Examples  include  counts  of  surface  imperfections  on  film,  flaws  in  fabric 
weave  and  nonconformities  in  completed  units  and  subassemblies.  The  par¬ 
ticular  application  developed  in  this  paper  relates  to  maintenance  activ¬ 
ities.  Frequently  maintenance  shops  process  similar  types  of  units,  such 
as  hydraulic  assemblies,  but  the  units  vary  substantially  in  size  and 
time  required  to  process  them.  In  such  cases,  it  may  be  reasonable  to 
assume  that  the  act  of  committing  an  error  in  processing  (the  >ccurence 
of  a  nonconformity)  has  a  constant  probability  as  a  function  of  t'ne. 

The  area  of  opportunity  for  the  occurrence  of  a  nonconformity  is  thus 
measured  in  units  of  time. 

We  assume  that  the  quality  control  procedure  consists  of  observing 
a  process  for  a  length  of  time,  H,  and  counting  the  number  of  noncon- 
formancts,  X,  that  occur  during  this  time  interval.  We  assume  that  X 


has  a  Poisson  distribution  with  intensity  A.  That  is,  the  mean  number 
of  nonconformances  observed  in  time  I!  is  AH.  Formulating  the  problem 
requires  the  specification  of  two  pairs  of  values: 

(1)  An  Acceptable  Process  Level,  A^,  and  its  associated 

risk  level  a.  A  Process  quality  level  that  is  con- 

*  *  o 

sidered  acceptable  as  a  process  average  measured  in  terms  of 
nonconformances  per  100  worker-hours.  The  probability  of  ac¬ 
cepting  the  hypothesis  that  the  process  is  operating  at  or  be¬ 
low  A  ,  when  it  actually  is  operating  at  A  ,  is  1  -  a. 
o  o 

(2)  A  Rejectable  Process  Level,  A ^ ,  and  its  associated 
risk  level,  6.  A^  is  the  process  quality  level  that  is  con¬ 
sidered  unacceptable.  The  risk  of  accepting  the  hypothesis 
that  the  process  is  operating  at  or  below  A^  when  it  actually 
is  operating  at  or  above  A^  is  6. 

The  two  points  (Aq,  1  -  a)  and  (A^g)  thus  define  the  operating  char¬ 
acteristic  curve  of  the  acceptance  control  chart  plan.  From  these  two 
points  we  will  derive  Che  Acceptance  Control  Limit,  K,  and  the  optimal 
subgroup  size,  H. 

Generally  speaking,  the  quality  control  procedure  will  involve  looking 

at  a  series  of  time  intervals,  H^,  ....,  and  observing  X^,  X^,  .... 

In  this  case,  we  assume  has  a  Poisson  distribution  with  mean  AH^.. 

til 

can  be  thought  of  as  the  size  of  t he  i  1  subgroup. 

The  intensity  of  nonconformances  at  the  acceptable  process  level 
(APL)  will  be  deonted  by  AQ,  while  the  intensity  at  the  rejectable  pro¬ 
cess  level  (RPl)  will  be  denoted  by  A^  With  K  denoting  the  acceptance 
control  limit  (ACL) ,  we  can  make  the  identifications  shown  in  Figure  1. 


Fig.  1  An  Acceptance  Control  Charting  Scheme  for  Poisson  Counts 

PROBLEM  SOLUTION 

Our  problem  is  to  determine  values  of  H  and  K  for  fixed  values  of  a, 

3*  and  A^.  Recall  that  we  want  to  choose  H  and  K  so  that  the  probabil¬ 
ity  of  X  exceeding  K,  when  A^  Is  the  true  intensity,  is  a;  and  the  proba¬ 
bility  of  X  being  less  than  or  equal  to  K,  when  A^  is  the  true  intensity,  is  3. 

We  shall  investigate  three  methods  of  calculating  H  and  K,  for  fixed 
Aq*  A^,  a  and  3.  The  first  will  use  exact  Poisson  probabilities;  the 
other  two  will  involve  normal  approximations  to  the  Poisson.  One  might 
ask  why  approximate  procedures  are  needed  when  an  exact  solution  is  known. 
The  answer  lies  in  the  fact  that  the  exact  solution,  for  all  possible  in¬ 
dustrial  applications,  requires  entensive  tables  of  Poisson  (or  Chi-square) 
probabilities.  Normal  approximations  have  been  used  in  industry  and  can 
be  worked  out  quickly  and  easily  with  reference  to  only  a  table  of  normal 


curve  areas. 


Exact  Solution 


We  can  find  the  exact  solution  for  H  and  E  based  on  Poisson  prob¬ 
abilities.  If  X  has  a  Poisson  distribution  with  mean  0,  then 
Pe(X  <  C)  .  P(X^(C  +  u  >  20) 

2 

where  denotes  a  Chi-square  random  variable  with  r_  degress  of  freedom. 
Thus,  a  table  of  Chi-square  probabilities  can  be  used  in  place  of  Poisson 
probabilities. 

Now  H  and  K  are  found  by  simultaneously  solving  the  equations 

and 

\*<*  i  «  ‘  p<x2«  +  1)  >  “iH)  '  B- 

These  equations  must  be  solved  iteratively. 


Standard  Normal  Approximation 

If  X  has  a  Poisson  distribution  with  mean  XH,  then  it  i9  well-known 

that 

X  -  XH 
/XH 


has,  approximately,  a  standard  normal  distribution  if  XH  is  large. 
(That  is,  (X  -  XH)//XH  has  a  distribution  which  tends  to  the  standard 
normal  distribution  as  XH  tends  to  infinity;  the  approximation  seems 
to  work  well  for  XH  greater  than  5.) 

If  z_  denotes  the  value  that  cuts  off  an  upper  tail  area  of  Y 
under  the  standard  normal  curve,  then  H  and  K  can  be  found  by  solving 
the  equations 


and 


z  = 
a 


K-XqH 


-z„ 

li 


K-X1H 

/TJF 
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Solving  these  equations  yields 


„  z  /X  +  z  /X, 
H  "  a  o  8  1 

X,  -  X 
1  o 


K  “  X.H  +  z  /XJH  -  X.H  -  z./Th 
0  a  0  1  6  1 


Square-Root  Normal  Transformation 

Since  X,  suitably  standardized,  is  approximately  normally  dis¬ 
tributed,  it  can  be  shown  that  /x  also  is  approximately  normally  dis¬ 
tributed.  The  variance  of  /x  is  essentially  free  of  X.»  for  large  aH, 
and  the  distribution  of  /X  tends  to  be  more  accurately  approximated 
by  a  normal  distribution  than  does  the  distribution  of  X,  for  moder¬ 
ate  values  of  XH. 

The  theory  (see  Johnson  and  Kotz  [1969])  actually  states  the 

2(/x  -  Ah) 


is  approximately  distributed  as  a  standard  normal  random  variable  if 
XH  is  large.  Working  on  the  true  square-root  scale,  we  find  H  and 
/IT  by  solving  the  equations 


z  =  2(/k~  -  /AaH  ) 
a  0 

and 

-zB  =  2(/F  -  /x^iT ) 

which  yield 

0.25  (z  +  z.)2 
„ - 2 - 1 

</*7-  /v 

and 

/lT  =  /XaH  +  1/2  z 
0  o 
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Transforming  back  to  the  original  count  scale 

K  -  (/ AnH  +  1/2  z  )2 
0  a 

The  equation  for  H  is  the  same  as  the  one  in  the  standard  normal  case 
when  a=g. 

NUMERICAL  STUDY 

The  values  of  H  and  K  were  found  for  various  fixed  values  of 
A  ,  A^,  o  and  6.  A  representative  sample  illustrative  of  our  findings  is 
shown  in  Table  1.  The  complete  set  of  results  is  contained  in  the  Appendix. 
To  illustrate  these  results,  we  will  look  at  the  first  row  of  figures 

where  a  =  0.025,  6  *  0.01,  Aq  =  0.1,  and  A^  =  0.6.  We  say  that  the 
process  is  out  of  control  if  there  are  more  than  K  nonconformances 
in  H  time  units  of  observation.  K  has  the  value  of  5  for  both  the 
standard  normal  approximation  and  the  exact  case  and  the  value  6  for 
the  square-root  normal  case.  Note  that  the  standard  normal  approxima¬ 
tion  gives  a  value  of  H,  23.469  time  units,  much  larger  than  the  true 
value,  22.000.  Thus  we  would  be  observing  the  process  longer  than 
we  should,  for  the  same  K  value,  and,  as  a  result,  have  a  greater 
probability  of  seeing  more  than  K  nonconformities  than  t he  nominal  value 
of  a  indicates. 

The  square-root  transformation  results  in  an  H  approximately  equal 
to  the  true  value,  but  the  K  is  slightly  larger.  Thus  the  probability 
of  seeing  more  than  K  nonconformities  would  be  slightly  smaller  than  the 
nominal  value.  This  pattern  prevails  throughout  most  of  the  cases  studied. 
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TABLE  1 


Values  of  H  and  K  for  Specified 
Aq»  X^,  a  and  B» 

a  -  0.025  6  *  0.01 


STANDARD 

NORMAL 

SQUARE  ROOT 

NORMAL 

FVAfT 

(CHI-SQUARE) 

Ao 

xi 

H 

K 

H 

K 

H 

K 

0.1 

0.6 

23.469 

5 

21.858 

6 

22.000 

5 

0.2 

0.7 

31.887 

11 

30.280 

12 

31.000 

11 

0.3 

0.8 

39.813 

18 

38.222 

19 

38.250 

18 

0.4 

0.9 

47.534 

27 

45.950 

28 

46.500 

27 

1.0 

4.0 

4.860 

9 

4.595 

10 

4.795 

9 

2.0 

5.0 

7.067 

21 

6.803 

22 

6.900 

21 

3.0 

6.0 

9.191 

37 

8.927 

38 

8.967 

37 

4.0 

7.0 

11.283 

58 

11.019 

58 

11.172 

58 

a  -  0.05 

B  *  0.01 

0.1 

0.6 

21.578 

4 

18.773 

5 

19.700 

4 

0.2 

0.7 

9 

10 

27.250 

9 

0.3 

0.8 

16 

16 

33.500 

15 

0.4 

0.9 

23 

39.442 

23 

39.556 

22 

1.0 

4.0 

4.408 

7 

3.944 

8 

4.000 

7 

2.0 

5.0 

6.299 

18 

5.839 

18 

5.860 

17 

3.0 

6.0 

8.131 

32 

7.663 

32 

7.767 

31 

4.0 

7.0 

9.915 

50 

9.459 

49 

9.536 

48 
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EXAMPLE  APPLICATIONS 


Example  1.  Check  of  Repaired  Items  Against  a  Standard.  The  data  of 

Table  2  shows  the  number  of  maintenance  errors,  X^,  observed  upon 

sampling  repaired  aircraft  parts  for  which  the  actual  repair  time  was 

hours.  In  this  example,  is  fixed  by  the  practical  sampling 

circumstances,  and  so  no  specific  A^  needs  to  be  determined.  It  is 

desired  that  A.  be  0.01  and  a  be  0.01.  Thus  z  is  2.33. 

0  -a 

TABLE  2  NONCONFORMANCES  AMONG  REPAIRED  ITEMS 


SAMPLE 

h 

Hi 

STANDARD  NORMAL 

K 

SQUARE  ROOT  NORMAL 
K 

1 

1 

58.33 

2.36 

3.72 

2 

4 

80.22 

2.89 

4.24 

3 

1 

209.24 

5.46 

6.82 

4 

2 

164.70 

4.64 

5.99 

Table  2  also  shows  the  values  of  K  obtained  by  the  standard  normal 
approximation  and  the  square-root  normal  transformation.  For  the  first 
sample : 

Standard  normal  approximation 

K  =  A„H  +  z  /TX 
0  a  0 

=  0.01(58.33)  +  2.33  /0. 01(58.33) 

=  2.36 

Square-root  normal  transformation 

K  =  (/XJT  +  1/2  za)2 
=  (/0. 01(58. 33)  +  2 . 33/ 2 ) 2 
=  3.72 
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Samples  1,  3  and  4  would  be  declared  "in  control"  at  the  standard 
value  of  Aq. under  either  scheme.  However,  sample  2,  with  =  4,  would 
be  declared  "out-of-control"  under  the  standard  normal  scheme  and  "in 
control"  if  the  square-root  normal  transformation  were  used;  the  observed 
value  is  very  close  to  the  boundary  in  either  case.  Whether  we  declare 

the  process  to  be  "out-of-control"  or  "in  control"  at  the  point  that 
sample  2  was  taken  depends  upon  whether  we  want  to  think  of  the  true  a 
risk  value  as  being  slightly  larger  than  0.01  or  slightly  smaller  than 
0.01.  In  many  cases  declaring  a  process  to  be  out-of-control  when,  in 
fact,  it  is  in  control  is  a  costly  error.  Thus  a  quality  control 
engineer  may  wish  to  use  the  more  conservative  procedure  that  lends 
itself  to  smaller  a  value. 

Example  2.  Establishing  a  Standard  Plan  to  Check  Maintenance  Errors 
in  a  Paint  Shop.  It  is  desired  to  set  up  a  standard  Acceptance  Control 
Chart  plan  for  checking  maintenance  errors  in  an  aircraft  subassembly 
paint  shop.  The  acceptable  process  level  is  3  errors  per  100  worker- 
hours  with  a  risk  level  (a)  of  0.05.  The  rejectable  process  level  is 
to  be  15  errors  per  100  worker-hours  with  a  risk  Level  of  0.10.  Values 
of  K  and  H  will  be  found  by  the  three  methods. 

As  previously  stated,  the  Chi-square  may  be  used  to  solve  for 
Poisson  probabilities.  Using  t he  Hald  Statistical  Tables  (1952); 

P<X2(K«)  i  2V»  '  “ 

P(X2(KH)  i  2‘lH>  *  1  -  6 

Substituting  the  values  of  A^,  A^,  a,  and  g  into  these  equations 
P(X2(K+1)  -  2<°-03)H)  =  °-05 
P(X2(K+1)  -  2(0-15),1)  =  U'9U 
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A  convenient  search  procedure  for  solving  these  equations  for  K  and  H 

2 

is  to  take  the  ratios  of  the  values  of  for  even  values  of  r  and 
solve  for  the  value  of  _r  that  is  closest  to  this  ratio.  The  value  of 

2 

H  may  then  be  found  from  the  resulting  values  of  Xf  taken  from  the 
table. 

2\ 

0-15/0.03-  5.0 

From  the  Hald  Tables  of  the  Chi-square  distribution: 

r  X2(r,  0.90)/X2(r,  0.05) 

6  10.6  /1.64  -  6.46 

8  13.4  /2 .73  =  4.91 

10  16.0  /3. 94  =  4.06 

Clearly,  the  ratio  of  the  two  Chi-squares  is  closest  to  the  desired 

value  of  5.0  when  r  equals  8.  The  value  of  K  then  must  be 
K  =  (8/2)  -1=3 

H  is  found  by  solving  the  equation 

2AH  =  X2 

r,  Y 

for  each  (A,  y)  pair  and  selecting  the  larger  (more  conservative  value.  Thus 

»  -  X(8,  0.90)  1  2\  -  UA  1  2(0'15)  ■  W'6'/ 

or 

H=X(8,  0i05)/2X0=  2.73/  2(0.03)  =45.50 

Thus  our  observation  time  should  be  45.50  hours. 

In  comparison,  using  the  standard  normal  approximation  yields  values 
of  H  and  K  of 

1 . 64  5  /  0~  1 6  3  +  1.282/0,15 
0,15  -  0,03 

K  =  0.03(42.41)  +  1.645/0. 03(42. 4l)~  =  3.13 
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By  the  square-root  normal  transformation,  these  values  are 

„  0.25(1.645  +  1.282)2  ,,  „„ 

H  =  - —  =  45. /3 

(/o.i5  -  /oToT) 

K  =  (/0. 03(46. 73)  +  1.645/2)2  -  4.03 

It  should  be  noted  that  the  actual  values  of  a  and  g  in  this  case 
are  0.040  and  0.122  using  the  standard  normal  approximation  and  0.014 
and  0.172  using  the  square-root  transformation.  Thus  both  approximations 
are  more  conservative  with  respect  to  a  error  and  less  conservative  with 
respect  to  8  error  than  the  plan  design  called  for  (a  =  0.05  and  8  =  0.10). 

As  with  many  cases  involving  observations  on  maintenance  operations 
the  actual  total  maintenance  time  involved  in  a  sample  subgroup  is  likely 
to  differ  from  the  planned,  or  design,  time.  Table  3  shows  the  actual 
subgroup  times  and  nonconformities  found  in  21  subgroups.  The  actual 
times  range  from  a  low  of  32.1  hours  to  a  high  of  57.8.  This  results 
from  the  fairly  wide  discrete  time  variation  required  to  process  a  unit. 

As  a  consequence,  it  may  be  necessary  to  recompute  the  control  limit 
based  upon  the  actual  time  in  a  subgroup  as  opposed  to  the  value  found 
for  the  design  time.  Since  the  count  of  nonconformities  is  integer¬ 
valued,  recalculation  of  the  control  limit  is  not  always  required. 

Figure  2  shows  the  Acceptance  Control  Chart,  using  three  sets  of 
control  limits,  for  the  sampling  data  of  Table  3.  Since  the  actual 
sample  hours  vary  from  subgroup  to  subgrcip,  it  is  inappropriate  to 
plot  a  central  line  on  this  chart.  (Where  the  sample  hours  can  be  held 
constant,  a  central  line  would  be  plotted  at  A^H.)  The  control  limit 
using  the  exact  Poisson  is  plotted  as  a  dash  line  at  the  value  3.5  for 
all  points  except  subgroups  7,  9,  12,  and  21.  Recalculation  was  necessary 
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32.1 


for  subgroups  7  and  12  (4.5)  because  of  the  larger  than  standard  sample 
hours  and  for  subgroups  9  and  21  because  of  the  smaller  than  standard 
sample  hours. 


ACL  by  exact  Poisson 

ACL  std.  normal  approximation 


Fig.  2  Acceptance  Control  Chart  for  Data  of  Table  3 

Table  3  also  shows  the  values  of  K  calculated  from  t he  standard 
normal  approximation  and  from  t he  square-root  normal  transformation. 
Again,  since  t he  count  of  nonconformities,  c,  is  integer-valued,  the 
ACL  for  each  approximation  has  been  set  half-way  between  integer ized 
values  of  K  and  K  +  l.  Where  these  ACL  values  differ  from  those  found 
by  the  exact  Poisson,  they  are  plotted  on  Figure  2  as  dotted  lines  for 
the  standard  normal  approximation  and  as  solid  lines  for  the  square-root 
normal  transformation. 

Observing  Figure  2,  it  should  be  noted  that  an  out-of-control 
condition  is  signalled  for  subgroup  1  by  both  the  exact  Poisson  and 
the  standard  normal  approximation  but  not  by  the  square-root  normal 
transformation.  In  those'  instances  wherein  the  ACL  by  the  standard 


normal  approximation  differs  from  the  exact  Poisson,  it  tends  to  be 


s?^Tf^ 


tighter.  Thus  the  standard  normal  approximation  tends  to  protect  more 
against  8  error  at  the  sacrifice  of  a  error.  The  square-root  normal 
transformation  tends  to  act  just  the  opposite.  Where  it  differs  from 
the  exact  Poisson  it  tends  to  be  looser  affording  greater  protection 
against  a  error  at  the  expense  of  8  error. 

This  feature  is  born  out  by  examination  of  Table  3  in  which  are 
tabulated  the  actual  a  and  8  error  for  each  subgroup  using  each  approxi¬ 
mation  method.  Recall  that  the  design  level  for  a  was  0.05.  By  the 
standard  normal  approximation,  actual  a  protection  ranged  from  0.035 
to  0.115  with  18  of  21  case  above  0.05.  In  the  case  of  the  square-root 
normal  transformation,  actual  a  error  ranged  from  0.015  to  0.039;  all 
cases  were  below  the  design  level  of  0.05. 

The  design  level  for  8  was  0.10.  In  the  case  of  the  standard  normal 
approximation,  the  actual  8  error  ranged  from  0.027  to  0.146  with  5  of 
21  cases  above  0.10,  For  the  square-root  transformation  all  but  three-v 
cases  were  above  the  design  level  with  the  actual  values  ranging  from 
0.067  to  0.292.  In  four  cases  the  actual  risk  levels  were  more  than 
double  the  design  level. 

It  should  be  noted  that  where  the  ACL  found  by  an  approximation 
method  agrees  with  that  found  by  the  exact  Poisson,  the  true  values  of 
a  and  8  apply  to  the  exact  Poisson  as  well.  Thus  when  actual  sample 
hours  differ  from  the  value  of  U  found  from  applying  the  Chi-square 
formulas,  the  actual  levels  of  protection  may  change  significantly. 

CONCLUSIONS 

This  paper  has  described  an  Acceptance  Control  Charting  approach 
for  process  control  of  cases  involving  the  observation  of  Poisson  counts. 
In  addition  to  deliniating  a  procedure  utilizing  the  exact  Poisson, 


procedures  for  use  of  the  popular  standard  normal  approximation  and 
of  the  not-so-frequently  used  square-root  normal  transformation  were 
developed  and  evaluated. 

It  was  shown  that  the  standard  normal  approximation  tended  to 
favor  protection  against  3  error.  To  the  extent  that  the  results  of 
ACL  calculations  differed  from  the  exact  Poisson,  the  difference  was 
biased  in  favor  of  3  error  protection.  On  the  other  hand,  usage  of 

the  square-root  normal  transformation  leads  to  ACL  calculations  offering 
better  protection  against  a  error.  To  the  extent  that  these  calculations 
differed  from  the  exact  Poisson,  the  bias  favored  a  error  protection. 

Study  of  a  number  of  cases,  of  which  Table  1  includes  a  sample, 
indicated  that  the  square-root  normal  transformation  gives  values  of 
H  and  K  that  oscillate  around  the  true  values  but  that  large  discrepancies 
between  the  approximate  and  true  values  are  rare.  We  therefore  recommend 
using  . ae  square-root  transformation  when  it  is  cumbersome  or  impossible 
to  use  exact  values  and  when  the  cost  of  a  error  is  high  in  relation  to 
3  error.  However,  for  those  cases  in  which  the  cost  of  3  error  is  equal 
to  or  greater  than  that  of  a  cror,  the  standard  normal  approximation 
is  preferable. 
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